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ABSTRACT 


Analysis  of  nine  events  recorded  by  the  NORSAR  short-period  array  reveals 
that  a beam  of  the  ten  Northeast  subarrays  has  only  .06  magnitude  units  less 
detection  capability  than  the  full  array  and  that  the  3C  subarray  has  only 
.18  magnitude  units  less  capability.  We  show  that  if  the  3C  subarray  were 
expanded  to  56  elements  inside  a 15  km  diameter  circle  it  would  nave  .3 
magnitude  units  more  detection  capability  than  the  present  full  array. 
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INTRODUCTION 


NORSAR  is  a combined  short-period,  long-period  22-subarray  seismic 
array  in  Norway,  depicted  in  some  detail  in  Figure  1.  Each  subarray 
contains  6 short-period  vertical  instruments  and  one  set  of  co-located 
Vertical,  North,  and  East  long-period  instruments  distributed  within  a 
10  km  diameter  circle.  Comprehensive  evaluations  of  the  short- 
period  NORSAR  array  have  been  published  by  Felix,  Gilbert,  and  Wheeler  (1971), 
Barnard  and  Whitelaw  (1972)  and  by  Ringdal  and  Whitelav  (1973).  These 
studies  have  covered  such  topics  as  variation  of  noise  spectrum  and  rms  noise 
l-vel  as  a function  of  time,  spectral  content  of  noise  and  signal,  optimal 
filtering,  signal  variation  across  the  array;  signal  loss,  noise  reduction, 
and  signal-to-noise  gain  in  beamforming,  detection  threshold  estimation, 
and  performance  of  short-period  discriminants  at  NORSAR. 

The  present  study  concentrates  on  the  question:  What  subarrays  and 

combinations  of  small  numbers  of  subarrays  will  give  a detection  threshold 
as  close  as  possible  to  that  of  the  full  array? 
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Figure  1.  Map  of  Norsar. 
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In  Table  I we  list  the  nine  events  considered  in  this  study.  They 
have  been  chosen  with  attention  to  their  distribution  in  distance  and 
azimuth,  and  source  regions  of  special  interest  have  been  included.  In 
addition  the  event  magnitudes  have  been  choser  so  that  there  is  no  clipping 
at  the  subarrays  with  the  largest  signal  amplitudes.  The  event  names  are 
the  same  as  those  used  by  Barnard  and  Whitelaw  (1972)  and  Ringdal  and 
Whitelaw  (1973). 

In  Figure  2 we  see  the  3C  and  12C  subarray  beams  for  KAZ/145/04N.  Each 
trace  is  shown  unfiltered,  filtered  0. 4-3.0  Hz,  and  filtered  1.1-2. 9 Hz.  The 
two  filters  are  3-pole  Butterworth  filters  created  by  program  FILCOF. 

We  see  from  subarray  12C  beams  that  the  1. 1-2.9  Hz  filter  improves 
the  signal-to-noise  ratio,  and  from  subarray  3C  that  the  signal  loss  is 
only  about  1.4  dB.  Investigations  by  all  of  the  authors  cited  in  the 
Introduction  to  this  report  reveal  that  a filter  pass  band  of  about  1. 0-3.0  Hz 
is  optimum  for  detection.  The  present  on-line  filter  for  beamforming 
NORSAR  data  is  a 3-pole  Butterworth  filter,  1.2-3. 2 Hz.  For  incoherent 
beamforming  the  on-line  filter  is  1.6-3. 2 Hz  (Tveitane,  1973).  We  shall 
perform  the  analysis  in  this  report  solely  with  the  1. 1-2.9  Hz  filter. 
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Figure  2.  Si.ourray  beams  for  KAZ/145/04N  for  subarravs  3C  and  12C  unfiltered, 
and  filtered  0. 4-3.0  and  1.1-2. 9 Hz.  Subarray  3C  shows  that  there  is  very 
little  signal  loss,  and  subarray  12C  shows  the  substantial  gain  in  S/N. 


SIGNAL  AND  NOISE  PROPERTIES  OF  THE  SUBARRAYS 


Using  program  QUICSAN  we  filtered  the  output  from  each  sensor  and  beamed 
each  subarray  using  the  azimuth  and  Herrin  (1968)  velocity  appropriate  for 
tne  event.  In  Table  II  we  see  the  results  for  each  subarray  for  each  event. 

The  signal  values  are  the  logarithm  to  the  base  10  of  half  the  maximum  peak- 
to-peak  subarray  beam  signal  in  millimicrons  in  Lhe  first  20  seconds  of  the 
signal.  The  noise  values  are  the  logarithm  of  the  root-mean-square  (rms) 
subarray  beam  noise  in  the  60  seconds  preceding  the  signal.  The  signal-to- 
noise  ratio  (S/N)  is  the  difference  of  these  two. 

In  Figure  3 we  see  these  three  numbers  contoured  on  a map  of  NORSAR  for 
the  event  KAZ/145/04N.  Note  that  the  contour  for  large  signal  and  the  contour 
for  small  noise  overlap  for  subarray  3C  which  has  the  largest  S/N  ratio. 

In  Figure  4 we  contoured  the  average  logarithm  of  the  noise  on  NORSAR  beams  for 
all  nine  events.  We  see  that  there  is  a maximum  noise  level  running  northwest- 
southeast  through  the  center  of  the  array;  and  by  reference  to  Figure  1 we 
can  see  that  it  is  associated  with  large  bodies  of  water  which,  as  might  be 
expected,  are  correlated  with  high  population  density  and  cultural  activity. 

The  noisiest  subarray,  14C,  can  be  seen  on  a more  detailed  map  to  lie  between 
two  roads  along  which  is  the  highest  density  of  population  on  the  map  outside 
of  a major  town.  Thus  it  would  appear  that  the  noise  levels  are  cultural  and 
can  amount  to  as  much  as  0.1-0. 2 magnitude  units. 

In  Figure  5 we  have  contoured  the  average  logarithm  of  the  signals 
on  NORSAR  beams  of  all  nine  events.  Arrows  indicate  the  direction  of  arrival 
of  the  nine  events.  (Events  from  the  Southwest  come  either  from  the  mid 
Atlantic  ridge  or  South  America.)  We  see  that  subarray  3C  is  the  best  subarray; 
and  that  the  northeast  corner  in  general  has  high  signal  levels. 

In  Figure  6 we  have  contoured  the  "difference"  of  Figures  4 and  5 
resulting  in  contours  of  S/N.  Again  we  see  that  the  northeast  subarrcys  are 
superior. 

In  Figure  7 we  have  plotted  the  name  of  the  event  next  to  the  subarray 
which  had  the  highest  S/N  for  that  event.  We  see  that  the  northeast  sub- 
arrays are  always  the  best,  except  for  the  RYU  event,  in  which  case  subarray 
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Figure  3.  Signal,  noise,  and  S/N  for  KAZ/145/04N.  Circle  at 
subarray  3C  indicates  size  of  the  subarrays. 


Log  signal  (0-P,  mp)  averaged  over  events  in  Table  I at  each 
Arrows  indicate  direction  of  arrival  of  events. 


Figure  5 
subarray 
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Figure  7.  Subarray  with  maximum  S/N  for  the  events  in  Table  I. 
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13C  is  best.  However,  it  is  only  0.04  magnitude  units  better  than  subarray 
4C.  This  shows  explicitly  that  the  superiority  of  the  northeast  subarrays 
is  not  due  to  a very  large  superiority  on  a few  events,  but  that  it  is  a 
general  phenomenon. 

The  general  results  found  in  Figure  5 are  similar  to  those  found  by 
Felix,  Gilbert,  and  Wheeler  (1971)  and  by  LaCoss  and  Filson  (1972).  Inspec- 
tion of  Figures  8 and  9,  which  have  been  drawn  from  data  published  by  these 
authors,  verifies  that  the  highest  signal  levels  are  found  in  the  northeast 
quadrant. 

Figure  10  (from  Berteussen,  1974)  shows  that  the  variability  in  travel 
time  residuals  as  a function  of  azimuth  is  a minimum  for  subarrays  2C  and  3C. 
This  suggests  that  the  geologic  structures  are  smooth  and  regular  under  these 
subarrays.  This  may  offer  some  explanation  for  their  high  signal  amplitudes. 
Berteussen,  Ringdal  and  Whitelaw  (1973)  attempted  to  explain  respectively  the 
travel  time  residual  and  amplitude  data  on  the  basis  of  an  irregular 
Mohorovicic  discontinuity.  All  solutions  in  Berteussen  (1974)  show  a 
relatively  flat  interface  under  2C  and  3C.  However,  crustal  measurements  by 
Kanestron  and  Haugland  (1971)  as  reported  by  Ringdal  and  Whitelaw  show  a 
sloping  interface  under  2C  and  3C. 


f 


Figure  10.  Length  of  arrow  is  proportional  to  B in  the  equation 
Relative  Residual  = A + Bsin0  where  0 is  the  azimuth  to  the 
event.  From  Berteussen  (1974). 
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SIGNAL  TO  NOISE  RATIOS  FROM  FULL  AND  PARTIAL  ARRAY  BEAMS 

To  beamform  the  full  array  it  is  necessary  to  use  travel  times  that  are 
corrected  to  the  world-wide  average  tables.  While  it  would  be  possible  to 
align  the  subarray  beams  precisely  for  these  large  events,  such  a 
procedure  would  not  be  a test  of  the  array's  operational  capability.  In 
routine  operations  it  is  necessary  to  use  a fixed  set  of  residuals,  and  if 
the  signal-to-noise  improvement  obtained  with  that  set  is  not  sufficient 
for  detection,  then  the  event  is  not  detected.  Any  further  improvement  which 
would  in  principle  be  possible  for  that  event  is  almost  academic.  An 
exception  to  this  would  be  if  an  event  is  detected  by  other  means  at  a 
site  very  close  to  a calibration  event.  Then  the  precise  residuals  for  that 
calibration  event  might  be  used  for  the  new  event.  In  general  it  will  be 
necessary  to  use  tabulated  corrections  for  routine  operations. 

The  most  comprehensive  table  of  such  corrections  for  NGRSAR  is  that 
published  by  Berteussen  (1974);  they  may  also  be  found  in  Anonymous  (1973). 
This  table  is  constructed  as  follows:  A set  of  one  of  more  calibration 

events  is  found  for  each  of  104  nodes  in  velocity  space.  The  known  geo- 
graphic locations  of  these  calibration  events  correspond  via  a known  world- 
wide travel  time  table  to  a location  in  velocity  space  (UCX,  UCY).  Each 
calibration  event  will  also  have  an  observed  location  in  velocity  space 
(UX,  UY)  which  is  obtained  for  each  calibration  event  by  fitting  a plane 
wave  to  the  arrival  times.  The  residuals  from  this  plane  wave  are  referred 
to  as  the  regional  corrections. 

To  determine  the  appropriate  delays  for  the  events  in  Table  I we  first 
found  the  closest  node  in  (UCX,  UXY)  space.  Then  the  vector  difference 
between  the  location  of  the  event  in  (UCX,  UCY)  space  and  the  location  of 
the  node  was  applied  as  a correction  to  the  location  of  the  node  in  (UX,  UY) 
space  to  give  the  location  of  the  event  in  (UX,  UY)  space.  Finally, 
we  used  the  regional  corrections  of  the  nearest  node.  Note  that  no  elevation 
corrections  are  required  since  they  are  absorbed  in  the  regional  corrections. 

The  delays  D(I)  are  computed  in  accordance  with  the  formulas  from 
Berteussen  (1974): 
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D(I)  = DPWF(I)  - DEV(I) 

DPWF(I)  = -(X(I)  • UX  + Y(l)  • UY). 


D(I)  is  the  delay  for  the  0 element  of  the  Ith  subarray,  DPWF  is  the 
moveout  appropriate  to  the  location  of  the  event  in  slowness  space,  X and  Y 
are  the  location  of  the  0 element,  and  DEV  are  the  regional  corrections.  In 
application  of  these  formulas  to  QUICSAN  it  is  important  to  remember  that 
all  time  delays  must  be  positive;  and  that  the  delay  for  a subarray  beam 
must  be  that  for  the  element  of  the  subarray  at  which  the  signal  first 
arrives.  The  delays  are  calculated  by  program  DELAY. 

In  Table  III  we  see  the  logarithms  of  the  signal,  noise,  and  signal-to- 
noise  for  the  full  array  beam;  a partial  array  beam  consisting  of  the  10 
subarrays  1-4B  and  1-6C;  and  subarray  3C  alone. 

We  note  that  the  average  signal  level  is  .15  magnitude  units  higher 
for  the  10-subarray  beam,  and  .43  magnitude  units  higher  for  subarray  3C. 

The  noise  level  is  .21  magnitude  units  higher,  and  .59  units  higher  respec- 
tively. The  net  result  is  a loss  in  S/N  of  .06  and  .18  magnitude  units 
respectively. 

This  last  result  implies  that  if  the  3C  subarray  were  expanded  to  14 
elements  its  detection  performance  would  equal  that  of  the  full  array. 
Carrying  this  one  step  further,  a 3C  subarray  of  56  elements  would  have  a 
detection  threshold  0.3  magnitude  units  lower  than  the  present  NORSAR.  With 
hexagonal  spacing  at  2 km  such  an  array  would  have  a diameter  of  15  km. 

To  analyze  the  question  of  whether  direct  summation  infinite  velocity 
subarray  beams  would  be  appropriate  for  NORSAR  we  evaluated  the  average 
loglO  (S/N)  for  infinite  velocity  beams  of  subarray  3C.  The  answer, 

1.56,  is  .45  m^  units  below  the  phased  subarray  beam  average  of  2.01 
seen  in  Table  III. 

In  this  analysis  we  included  the  two  close-in  events  from  the 
Urals  and  Greece,  URA  and  GRE,  at  distances  of  20  and  23  degrees.  As  would 
be  expected,  the  travel  time  residuals  and  signal  correlation  between 
subarrays  are  poor  for  these  events.  Thus  the  expanded  3C  array  should  be 
even  more  superior  for  these  events.  However,  excluding  these  events  results 
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in  rj  significant  change  in  the  averages.  Some  of  the  loss  of  detection 
capability  of  the  full  array  for  regional  events  is  presently  recovered  in 
practice  at  NORSAR  by  incoherent  beamforming;  see  for  example  Ringdal, 

Husebye  and  Dahle,  1972;  and  Blandford  and  Wirth,  1973. 

Table  IV  shows  that  the  average  noise  reduction  by  the  full  array  beam, 
and  by  the  10-subarray  beam,  with  respect  to  the  average  individual  sensor 
noise  level,  is  proportional  to  the  square  root  of  the  total  number  of 
sensors  in  the  beam.  The  signal  loss  for  teleseismic  events  with  reject 
to  individual  sensors  is  4.9  and  4.2  dB  respectively;  and  the  signal-to-noise 
ratio  gain  is  16.0  and  12.4  dB  respectively. 

These  results  on  signal-to-noise  ratio  improvement  could  conceivably  be 
misleading  if  the  signal-to-noise  ratio  were  distributed  log-normally,  as 
numerous  investigators  have  shown  to  be  the  case  at  LASA  and  NORSAR,  and  if 
the  variance  at  NORSAR  were  especially  large.  This  would  be  true,  because  if 
the  variance:  were  large,  the  upper  tail  of  the  log-normal  distribution  would 
heavily  weight  the  "average"  S/N  ratio.  The  average  would  then  be  much 
higher  than  the  value  of  the  median  or  typical  subarray  which  is  the  measure 
of  improvement  most  in  accord  with  the  question,  "What  is  the  improvement  of 
a large  array  over  a typical  small  array?  . 

To  examine  this  problem  we  constructed  Table  V which  shows  the  standard 
deviation  in  magnitude  units  of  signal,  noise  and  S/N  among  individual 
elements  for  the  10  events.  We  see  that  for  teleseismic  events  the  signal  loss 
relative  to  the  average  individual  sensors  is  4.9  dB  using  an  arithmetic 
average,  and  3.6  dB  using  a logarithmic  average.  The  difference,  1.3  dB, 
would  not  be  significant  in  most  practical  array  or  network  design  problems. 

As  a side  light  it  is  surprising  that  the  two  regional  events  have  the 
smallest  standard  deviation  for  the  signal,  even  though  we  saw  in  Table  IV 
that  the  array  beam  signal  loss  is  large  for  these  events  due  to  poor 
signal  correlation.  Perhaps  this  shows  that  the  regional  signals  are  so 
thoroughly  diffused  that  the  average  maximum  amplitudes  are  more  scab.e 
than  they  are  for  nore  vertically  arriving  rays  in  which  sharp  focur  ig  and 
defocusing  effects  are  predominant. 
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SUMMARY  AND  DISCUSSION 


• Noise  levels  vary  across  the  array  by  0.1-0. 2 magnitude  units  and  seem 
to  be  correlated  with  the  level  of  cultural  activity. 

• High  signal  levels  from  events  at  all  azimuths  are  recorded  by  the 
northeast  subarrays. 

• The  northeast  subarrays  have  high  S/N  ratios. 

• An  array  of  10  northeast  subarrays  has  only  .06  magnitude  units  less 
detection  capability  than  the  full  array. 

• The  3C  subarray  alone  has  only  0.18  magnitude  units  less  detection 
capability  than  the  full  array. 

• If  the  3C  subarray  were  expanded  to  a 56  element,  15  km  diameter  array 
with  2 km  spacing,  its  detection  threshold  would  be  0.3  magnitude  units 
lower  than  the  present  full  array.  Such  a small  array  would  require  much  less 
computer  power  for  analysis  than  does  the  full  array  and  would  better  preserve 
the  high-frequency  energy  which  may  be  useful  for  discrimination  purposes.  A 
possible  method  of  implementation  would  be  to  move  the  equipment  from  existing 
southwest  subarrays  into  sites  around  subarray  3C. 
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